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1. Introduction 

The search for New Physics (NP) together with stringent tests of the Standard Model (SM) 
pass through a detailed study of physical processes involving the b-quaik. Two particularly im- 
portant cases for the detection of potentially large NP contributions are the purely leptonic decays 
B — > x v T and B s — > pL~ /i + . The relevant entries in the SM prediction for these decay rates are 
the CKM matrix element V u b, and the pseudoscalar decay constants fs and fg s . Additionally the 
B-Physics parameters, Am^j Am s , Am^, play a crucial role in the Unitarity Triangle Analysis, and 
their determination rely on lattice QCD computations of the bag -parameters Bb s and BbJBbj- 

Despite the present (B-factories, LHCb) and future (SuperB factories) experimental programs, 
lattice QCD results of hadronic parameters need to have reduced uncertainties at the level of ~ 1 % 
(see refs [2, 3] for recent reviews). 

In this proceeding we report on the ongoing project, within the ETMC collaboration, to com- 
pute B-physics hadronic parameters with Wilson twisted mass (Wtm) lattice fermions. We extend 
and improve our previous analysis [4, 5] in two ways: 

• we optimize the interpolating operators for heavy-light systems to better project onto the 
fundamental state 

• we extend the range of heavy masses, /X/,, considered reaching values of approximately fi^ ~ 
2.5 m c , where m c is the charm quark mass. 

2. Computational details 

We use, for this analysis, the Nf = 2 dynamical gauge configurations with up and down mass de- 
generate quarks, m u u = pin generated by the European Twisted Mass Collaboration (ETMC) [6]. 
The lattice action is the tree-level improved Symanzik gauge action [7] and the twisted mass quark 
action [8] at maximal twist [9]. The strange and the charm quarks are quenched in this work. 
We have used four lattice spacings [10] and we take the values of the light and strange renormal- 
ized quark masses together with the pseudoscalar density renormalization constants Z^* s (2 GeV) 
respectively from [10] and [11]. All the lattice details of our computation can be found in our 
recent lattice 2012 contribution [12]. As usual heavy-light meson masses and decay constants are 
extracted studying the Euclidean time dependence of the 2-point functions with heavy-light interpo- 
lating operators. To improve the projection onto the fundamental state and keep the noise-to-signal 
ratio under control we extract meson masses at relatively small temporal separations. Thus the 
lowest lying energy eigenstate must have a large overlap with our interpolating fields. We have 
constructed these fields using the so-called Gaussian smearing [13] (see [12] for details). The im- 
provement with respect with the standard local interpolating operators can be seen in the left plot 
of fig. 1 where we show the Euclidean time dependence of the effective masses for a particular sim- 
ulation point obtained from local-local (LL), and two improved interpolating operators. The plot 
describes well the general properties for heavy-light effective masses: the improved correlators 
allow a safe extraction of the ground state mass at shorter Euclidean time separations compared to 
the LL correlators. 

'Very recently the first measurement of this decay rate [1] shows a good agreement with the SM prediction. 
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Figure 1: Left plot: Euclidean time dependence of the effective masses at a lattice spacing of a = 
0.098(3) fm, of a heavy-light correlation function obtained from LL and improved correlation functions. 
The heavy-quark mass corresponds to ~ 2.3m c . Right plot: chiral-continuum extrapolation of the ratio of 
heavy-light meson masses at the heaviest quark masses, jui ' . 



3. The £-quark mass and decay constants fg and fs s 

To determine the B-physics quantities we implement the ratio-method (see ref. [5] for the 
details of the method). We shortly recall here the basic idea for the computation of the ft-quark 
mass. The same strategy applies for the other quantities computed in these proceedings. The static 
limit of the heavy-light meson mass M^i in the pole heavy-quark mass nf° le is given by 

lim (m^/aC 16 ) = 1 . (3.1) 

We take now a set of heavy-quark masses juj 1 ^ < /l[ 2 ' < • • • < jur^ with ratio A, i.e. fir = Apr % 

fixed. We denote the masses renormalized in the MS at a 2 GeV scale with a "bar". The ratios of 

( n \ 

heavy-light meson masses at subsequent values of the heavy-quark mass, y(n h , A; Hi, a), properly 
normalized, are the key quantities of the ratio-method [5]. From QCD asymptotic freedom and 
eq. (3. 1) it follows that the ratio y in the continuum limit 

Jim lim^pj^, Ajp^c?) = 1 , (3.2) 

has an exact static limit. The chiral and continuum extrapolation of the ratio of the heavy-light 
meson masses evaluated at the heaviest quark mass is shown in the right plot of fig. 1. The lattice 
results are well described by a fit linear in a 2 and in the light quark mass. We also note that 
discretization errors are sufficiently well under control. As it can be seen from the left plot in 
fig. 2, the ratio y(p /n A), extrapolated to the chiral and continuum limit, has a non-perturbative 
heavy-quark mass dependence that is well described by the HQET-inspired function 

,(ft,A) = l + ^ + ?#. (3.3) 

This plot clearly shows that the accurate description of the non-perturbative dependence of the ratio 
v on Jl h is possible only because of the improvements that we have implemented in this analysis. 
In particular the results are well described by the fit ansatz in eq. (3.3), that provides, with the fit 
parameters 77 1 and T72, the non-perturbative description of the heavy-quark mass dependence of the 
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Figure 2: Left plot: the continuum ratio-function y(ju A , X ), for X — 1 . 1784 as a function of the heavy-quark 
mass. For comparison we show our previous results [4] without the interpolating operators improvement. 
Right plot: chiral extrapolation for the double-ratio [/)„/ //,/] ■ [/ s // /^ xp ] at the triggering point, ju^ 1 . 

ratio y(p !v X) over the whole range of heavy-masses. The Z?-quark mass can be computed simply 
taking the triggering heavy-light meson mass M#(p^ ) around the charm region and through a 
repeated application of the ratio y(JI h ,X) reaching the experimental value of the B meson mass M#. 
This is always possible with a slight tuning of X and that leads after a 4-loop evolution to 2 

mf"{m b ) \ Nf=2 =4.35(12) GeV. (3.4) 

We obtain a perfectly consistent result if we use as input the heavy-strange meson mass M\ ls . The 
result we present here (3.4) is in agreement with our previous determination [4] with slightly re- 
duced eiTor. The 2.7% relative error is dominated by the ones associated to the lattice spacing and 
Zp determinations. The statistical and systematic uncertainties coming from the application of the 
ratio method turn out to be negligible, after the improvements just discussed. 

A completely analogous strategy can be adopted to compute heavy-light and heavy-strange de- 
cay constants. We find advantageous to define ratios with exactly known static limit for the heavy- 
strange decay constant, z s (W^ ,X;Jli,Jl s ,a), and for the ratio /^//g, ^(Ji^\X\~\i h ~\i s ,d) [4]. Both 
fhs(Ph ) at tne triggering point and the ratio z s have a smooth chiral-continuum extrapolation with 
cutoff effects always well under control. The heavy-quark mass dependence in the continuum of z s 
can again be described by a formula as the one in eq. (3.3). We use the ft-quark mass, m^, we have 
determined as input leading us, using the ratio-method, to 

f Bs = 234(6) MeV. (3.5) 

The double-ratio £(jU/j , A;/X;,/X 5 ,a) has a very smooth and weak dependence on the light- 
quark mass, the lattice spacing and the heavy-quark mass. The chiral behaviour of fh s / fhi at the 
triggering mass, fr, % suggests that the heavy-quark still behaves as a relativistic one. Thus for 
the chiral-continuum extrapolation of fhs/fhl(P^) we find convenient to study the double-ratio 
[fhs/fhi] ' [fsi/ft? P ] > f° r which a smooth dependence on the light quark mass is expected [14, 15]. 
The results are well described by a linear fit in the light-quark mass as it can be seen from the right 

2 Note that a fc-quark mass value smaller by ~ 120 MeV would be found by setting Nf = 4 (rather than 2) in the 
evolution from 2 GeV to the fc-scale. 
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plot in fig. 2. This analysis of > A'.jll^jllJ and [fhs/fhi] • [/k//g Xp ] leads to the values 

^ = 1.19(5), / B = 197(10) MeV, (3.6) 
is 

where we add in quadrature the statistical and systematic uncertainties. 
4. Bag parameters 

The renormalized operator 0^ =1 = [by^ (1 — Ys)q] [by^ (1 — 75)9] in QCD is parametrized in 
terms of the bag parameter Bb ci (fi) as follows 

(S,|Of = 2 \B q ) m = ^f B Bl q ^)Ml q (4.1) 

where is the renormalization scale. With Wtm fermions one can define an operator with the cor- 
rect continuum limit whose renormalization is multiplicative [16]. Using heavy-quark symmetry 
one can show that the ratio of renormalized B-parameters evaluated in QCD is expected to approach 
unity as 1 /Jl h — > 0, with corrections of order 1 /\og(ji h / Aqcd) for small values of 1 /Jl h . To analyze 
the impact of such corrections we consider the ratio 

B Bv (Ji h ,JI h a;n) Cfe;M*,A0 

where the C-factors ratio contains the info on the 1 /log(/I /l )-cori - ections at a fixed order in renor- 
malization group (RG) improved perturbation theory (PT). Matching HQET to QCD in PT one 
can evaluate the function C. We consider here HQET-to-QCD matching only at tree level and 
LL order in PT (thereby avoiding the complications of operator mixing in HQET [17]), and con- 
firm the impact of 1 /log Ql/j) -corrections on the final results to be at the level of one standard 
deviation [18]. The chiral and continuum limit of the ratio (O s at the heaviest mass is well under 
control, see left plot of fig. 3. We refer to our contribution [18] to the lattice 2012 conference 
for additional details. The heavy-quark mass dependence of (O s is well described by the formula 
(O s (Ph) = 1 + c i W/M/i + c 2(A)//I^, see right plot of fig. 3. Applying the ratio-method one finally 
gets 

Bf{m h )= 0.90(5), (4.3) 

where the error is the sum in quadrature of statistical and systematic uncertainties. If we apply the 
ratio-method to Bs d and BbJBbj we obtain the preliminary estimate of 

= 1.03(2) , Bl s (m b ) = 0.87(5) . (4.4) 

BB d 

Given this result we can also give a preliminary estimate for the parameter 

e fs s JBb 



/bvBb 

where the error is the sum in quadrature of all uncertainties. 



1.21(6), (4.5) 



5 



B-physics from lattice QCD 



A. Shindler 




MGeV] 

"0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 



1.08 



1.06- 



1.04- 



1.02- 



0.98, 





1 


1 1 










1 


l/H h [GeV'] - 



0.6 



Figure 3: Left plot: chiral-continuum extrapolation of the ratio of the bag parameters, co s , for the heaviest 
quark mass analyzed, fj, h . Right plot: heavy-quark mass dependence of the ratio function 03 s (fl h ,X), for 
A = 1.1784 in the continuum and its non-perturbative description. 
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